Cellular Morphology and Transcriptome Comparative Analysis of Astragalus membranaceus Bunge Sprouts Cultured In Vitro under Different LED Light

Astragalus membranaceus, the major components of which are saponins, flavonoids, and polysaccharides, has been established to have excellent pharmacological activity. After ginseng, it is the second most used medicinal plant. To examine the utility of A. membranaceus as a sprout crop for plant factory cultivation, we sought to establish a functional substance control model by comparing the transcriptomes of sprouts grown in sterile, in vitro culture using LED light sources. Having sown the seeds of A. membranaceus, these were exposed to white LED light (continuous spectrum), red LED light (632 nm, 1.58 μmol/m2/s), or blue LED light (465 nm, 1.44 μmol/m2/s) and grown for 6 weeks; after which, the samples were collected for transcriptome analysis. Scanning electron microscopy analysis of cell morphology in plants exposed to the three light sources revealed that leaf cell size was largest in those plants exposed to red light, where the thickest stem was observed in plants exposed to white light. The total number of genes in A. membranaceus spouts determined via de novo assembly was 45,667. Analysis of differentially expressed genes revealed that for the comparisons of blue LED vs. red LED, blue LED vs. white LED, and red LED vs. white LED, the numbers of upregulated genes were 132, 148, and 144, respectively. Binding, DNA integration, transport, phosphorylation, DNA biosynthetic process, membrane, and plant-type secondary cell wall biogenesis were the most enriched in the comparative analysis of blue LED vs. red LED, whereas Binding, RNA-templated DNA biosynthetic process, DNA metabolic process, and DNA integration were the most enriched in the comparative analysis of blue vs. white LED, and DNA integration and resolution of meiotic recombination intermediates were the most enrichment in the comparison between red LED vs. white LED. The GO term associated with flavonoid biosynthesis, implying the functionality of A. membranaceus, was the flavonoid biosynthetic process, which was enriched in the white LED vs. red LED comparison. The findings of this study thus indicate that different LED light sources can differentially influence the transcriptome expression pattern of A. membranaceus sprouts, which can provide a basis for establishing a flavonoid biosynthesis regulation model and thus, the cultivation of high-functional Astragalus sprouts.


Introduction
Astragalus membranaceus Bunge, commonly referred to as Mongolian milkvetch, is an herbaceous perennial plant belonging to the Leguminosae family. In Republic of Korea, high-quality A. membranaceus is grown in Jeongseon, Gangwon-do, and has been used as a medicinal plant since ancient times [1,2]. A. membranaceus is naturally widely distributed in Asian regions such as Korea and China and in parts of Europe and Africa [3]. In oriental medicine, it is used for diuresis, tonicity, and blood pressure lowering and has also been reported to have blood sugar regulation, antitumor, and antiviral activities [4]. A representative component contained in A. membranaceus is formononetin, an isoflavone glycoside, which as a phytoestrogen is established to be a natural substitute for female hormones [5]. In addition, saponins, flavonoids, amino acids, trace elements, and polysaccharides have been reported as major biologically active constituents of A. membranaceus [6].
Consequently, environmental damage exacerbated by ongoing climate change and plant factories is increasingly gaining attention as an alternative approach for cultivating high-quality crops in a fixed environment without being influenced by the external environment [7]. Although compared with traditional field cultivation, there are certain disadvantages associated with factory-based cultivation, notably the higher equipment and maintenance costs; various studies are being conducted with the aim of establishing systems that compensate for these disadvantages [8]. Cultivation under plant factory conditions is based on the mechanical regulation of light, humidity, and temperature, and in this regard, there are an increasing number of studies focusing on the use of light-emitting diodes (LEDs) [9]. LEDs are mercury-free, environmentally friendly, and lightweight sources of light that are energy efficient, have a long lifespan, and have advantages such as simple circuitry and the provision of light of a specific quality [10]. LED light can be modified to produce light of different qualities using the principle that light is generated as electron flows, and the light quality can be selected according to the user's needs [11]. As such, the use of LEDs with specific light quality has been assessed with respect to the cultivation of a range of vegetable crops, although it has rarely been assessed for the cultivation of medicinal plants.
Although the basic data relating to the utilization of LEDs for plant cultivation is gradually accumulating, information pertaining to the identification and control of the relationships between LED light sources and the germination, growth, and contents of functional substances is still far from complete. To facilitate the bulk production of functional substances in crop plants based on the control of LED lighting, it will be necessary to investigate the influence of these light sources on the pattern of genes associated with functional component biosynthesis [12]. Based on this biosynthesis, it would be possible to modify different biological activities and thus, the contents of the functional substances of crops. In this regard, next-generation sequencing (NGS) technology can be applied to obtain large-capacity fragment sequencing data more rapidly and at a lower cost compared with traditional sequencing methods. NGS-based RNA sequencing (high-throughput mRNA sequencing, RNA-seq) can be used to sequence transcriptomic cDNA and quantify the relative amounts of transcript expression. It can be used to identify genes of interest that are specifically expressed according to research objectives, and offers the potential to analyze gene expression, even in the absence of relevant reference genome information [13].
Transcriptomic analysis is a particularly important approach for the comparison of biological activities and differences in gene expression in different plant tissues and for elucidating genome functionality, and numerous high-quality new crop varieties have been developed based on the genetic improvement of the plant's secondary metabolites using genetic engineering technology. Moreover, advances in transcriptomics have made it possible to identify genes associated with specific functions and predict hitherto unknown genes. Recently, the scope of this research expanded to include medicinal plants, for which comparative transcriptome analyses have been successfully performed [14,15]. To complement this growing body of research, we sought in the present study to analyze the transcriptomic responses of A. membranaceus sprouts grown under in vitro conditions in which plants were exposed to LED light sources of different wavelengths. This comparative transcriptome analysis enabled us to establish a light environmental control model that could be applied to enhance the production of functional substances of interest in sprouts of the medicinal plant A. membranaceus.

Comparison of the Cellular Morphologies of A. membranaceus Sprouts Germinated In Vitro under Three Different LED Light Sources
To compare the tissue morphologies of in vitro-cultured A. membranaceus sprouts exposed to different light sources for 6 weeks, we performed scanning electron microscopy observations of leaf and stem cross-sections ( Figure 1). Among the different treatments, the largest leaf cells (58.57 ± 6.17 µm) were observed in plants exposed to the red LED, whereas plants illuminated with white LED light were found to have the smallest leaf cells (20.67 ± 6.50 µm) by cross-sections. With regards to stem cell width, the highest and lowest values of 22.17 ± 2.44 and 16.81 ± 1.89 µm were observed in plants exposed to white and red LED lights, respectively (Table 1). [14,15]. To complement this growing body of research, we sought in the present study to analyze the transcriptomic responses of A. membranaceus sprouts grown under in vitro conditions in which plants were exposed to LED light sources of different wavelengths. This comparative transcriptome analysis enabled us to establish a light environmental control model that could be applied to enhance the production of functional substances of interest in sprouts of the medicinal plant A. membranaceus.

Comparison of the Cellular Morphologies of A. membranaceus Sprouts Germinated In Vitro under Three Different LED Light Sources
To compare the tissue morphologies of in vitro-cultured A. membranaceus sprouts exposed to different light sources for 6 weeks, we performed scanning electron microscopy observations of leaf and stem cross-sections ( Figure 1). Among the different treatments, the largest leaf cells (58.57 ± 6.17 μm) were observed in plants exposed to the red LED, whereas plants illuminated with white LED light were found to have the smallest leaf cells (20.67 ± 6.50 μm) by cross-sections. With regards to stem cell width, the highest and lowest values of 22.17 ± 2.44 and 16.81 ± 1.89 μm were observed in plants exposed to white and red LED lights, respectively (Table 1).     (Table 2). Based on de novo assembly, we detected a total of 45,667 genes in the germinated sprouts.  Table 3 shows the mapping results, which were applied to analyze the genes expressed in A. membranaceus exposed to the three different light sources, using the 45,667 genes identified via de novo assembly as a reference. For plants cultivated under blue, red, and white LED lights, we obtained 25,740,782 (70.1%), 26,850,570 (68.4%), and 20,898,430 (65.9%) mapped reads, respectively, with an average of 68.1% of the data being successfully mapped to the reference genome. Based on the confirmation of gene expression levels via differentially expressed gene (DEG) analysis, we established the numbers of genes to which at least one clean read was mapped for each of the 45,667 total genes, with 42,572; 42,542; and 42,606 clean data being mapped for plants exposed to blue, red, and white LED light, respectively (Table 4). Using the criterion of a log 2 -fold change, we established that the number of up-and downregulated genes between blue LED and red LED, blue LED and white LED, and red light and white LED treatment groups were 132 and 153, 148 and 93, and 144 and 91, respectively (Table 5). Figure 2 shows an MA Plot generated using the data presented in Table 5, in which genes showing significant differences in expression between two samples are indicated by a red color. To confirm correlations between the two samples, we generated a correlation plot, on the basis of which we determined that the highest correlation between different treatments (0.94) was obtained for A. membranaceus germinated in vitro under white and blue LED light ( Figure 2).

GO Analysis of A. membranaceus Sprouts Germinated In Vitro under Three Different Light Sources
The results of GO analysis of up and downregulated genes revealed a significant difference in expression level when comparing the two samples based on DEG analysis (Figures 3 and 4). When the same treatment group was compared, cellular process, catalytic activity, and cellular process showed the highest expression levels with 51, 50, and 56 genes among the upregulated genes, respectively. For blue LED vs. red LED and blue LED vs. white LED comparisons, the largest proportions of downregulated genes were found to be involved in catalytic activity (68 and 32 genes, respectively). For the comparison of red LED and white LED treatment groups, we found that 33 genes associated with cellular processes were the most downregulated genes. A comparison of the expression levels of up-and downregulated genes in the different GO functional categories revealed that in the molecular function category, genes associated with binding and catalytic activity were characterized by the highest levels of expression. In the cellular

GO Analysis of A. membranaceus Sprouts Germinated In Vitro under Three Different Light Sources
The results of GO analysis of up and downregulated genes revealed a significant difference in expression level when comparing the two samples based on DEG analysis (Figures 3 and 4). When the same treatment group was compared, cellular process, catalytic activity, and cellular process showed the highest expression levels with 51, 50, and 56 genes among the upregulated genes, respectively. For blue LED vs. red LED and blue LED vs. white LED comparisons, the largest proportions of downregulated genes were found to be involved in catalytic activity (68 and 32 genes, respectively). For the comparison of red LED and white LED treatment groups, we found that 33 genes associated with cellular processes were the most downregulated genes. A comparison of the expression levels of up-and downregulated genes in the different GO functional categories revealed that in the molecular function category, genes associated with binding and catalytic activity were characterized by the highest levels of expression. In the cellular component category, except for the downregulated genes between the blue LED and white LED treatment groups, we detected high levels of expression for all cellular anatomical entity-related genes. With respect to the biological process category, we detected generally high expression levels among genes associated with cellular processes. Expression levels of metabolic processrelated genes were also high, with the exception of the downregulated genes between blue LED and white LED treatment groups and upregulated genes between the red LED and white LED groups. component category, except for the downregulated genes between the blue LED and white LED treatment groups, we detected high levels of expression for all cellular anatomical entity-related genes. With respect to the biological process category, we detected generally high expression levels among genes associated with cellular processes. Expression levels of metabolic process-related genes were also high, with the exception of the downregulated genes between blue LED and white LED treatment groups and upregulated genes between the red LED and white LED groups.

qPCR Analysis for Reference Genes by GO Analysis
To confirm the reliability of the transcriptome comparative data, qPCR was performed by selecting five genes upregulated in white LED light and five genes upregulated in red LED light based on the blue LED light ( Figure 5). As a result of the qPCR analysis, it was confirmed that the expression patterns of all genes upregulated in the white and red light sources were represented by correlating in the transcriptome analysis. In particular, among the 10 genes, it was confirmed that the expression levels of the component category, except for the downregulated genes between the blue LED and white LED treatment groups, we detected high levels of expression for all cellular anatomical entity-related genes. With respect to the biological process category, we detected generally high expression levels among genes associated with cellular processes. Expression levels of metabolic process-related genes were also high, with the exception of the downregulated genes between blue LED and white LED treatment groups and upregulated genes between the red LED and white LED groups.

qPCR Analysis for Reference Genes by GO Analysis
To confirm the reliability of the transcriptome comparative data, qPCR was performed by selecting five genes upregulated in white LED light and five genes upregulated in red LED light based on the blue LED light ( Figure 5). As a result of the qPCR analysis, it was confirmed that the expression patterns of all genes upregulated in the white and red light sources were represented by correlating in the transcriptome analysis. In particular, among the 10 genes, it was confirmed that the expression levels of the

qPCR Analysis for Reference Genes by GO Analysis
To confirm the reliability of the transcriptome comparative data, qPCR was performed by selecting five genes upregulated in white LED light and five genes upregulated in red LED light based on the blue LED light ( Figure 5). As a result of the qPCR analysis, it was confirmed that the expression patterns of all genes upregulated in the white and red light sources were represented by correlating in the transcriptome analysis. In particular, among the 10 genes, it was confirmed that the expression levels of the reference genes Gene_174710T and Gene_334410T, which are upregulated in the white LED light source, were significantly higher than those of other LED light sources. reference genes Gene_174710T and Gene_334410T, which are upregulated in the white LED light source, were significantly higher than those of other LED light sources.

Up-and Downregulation of Transcriptomes in A. membranaceus Sprouts Treated with Blue LED light vs. Red LED light
The GO term with the highest number of upregulated transcriptomes was binding (GO:0005488), for which 12 were counted. The GO terms with the next highest number of upregulated transcriptomes (six) were DNA integration and transport. The GO terms with five counted transcriptomes were phosphorylation, DNA biosynthetic process, membrane, and plant-type secondary cell wall biogenesis. There were 12 transcriptomes counted as two and 28 transcriptomes counted as one (Table 6). Among the downregulated transcriptomes, the most represented GO term was binding (GO:0005506) with eight counts, followed by DNA integration (GO:0015074) with seven counts. The five-count GO term was the integral component of the membrane, and the four-count GO term was carbon utilization. The three counts were obtained for telomere maintenance (GO:0000723), defense response (GO:0006952), regulation of DNA-templated transcription (GO:0006355), protein ubiquitination (GO:0016567), and methylation (GO:0032259). Two counts were classified as GO terms for 11 transcriptomes, and one count was classified as the GO term for 50 transcriptomes (Table 7). Table 6. GO enrichment of upregulated transcripts in blue LED vs. red LED light of A. membranaceus sprouts (p-value < 0.05).

Up-and Downregulation of Transcriptomes in A. membranaceus Sprouts Treated with Blue LED Light vs. Red LED Light
The GO term with the highest number of upregulated transcriptomes was binding (GO:0005488), for which 12 were counted. The GO terms with the next highest number of upregulated transcriptomes (six) were DNA integration and transport. The GO terms with five counted transcriptomes were phosphorylation, DNA biosynthetic process, membrane, and plant-type secondary cell wall biogenesis. There were 12 transcriptomes counted as two and 28 transcriptomes counted as one (Table 6). Among the downregulated transcriptomes, the most represented GO term was binding (GO:0005506) with eight counts, followed by DNA integration (GO:0015074) with seven counts. The five-count GO term was the integral component of the membrane, and the four-count GO term was carbon utilization. The three counts were obtained for telomere maintenance (GO:0000723), defense response (GO:0006952), regulation of DNA-templated transcription (GO:0006355), protein ubiquitination (GO:0016567), and methylation (GO:0032259). Two counts were classified as GO terms for 11 transcriptomes, and one count was classified as the GO term for 50 transcriptomes (Table 7).

Up-and Downregulation of Transcriptomes in A. membranaceus Sprouts Treated with Blue LED Light vs. White LED Light
Among the transcriptome classified as upregulated, the most represented GO term was Binding (GO:0005488) with eight counts. RNA-templated DNA biosynthetic processes were classified as a six-count transcriptome, and five counts were obtained for DNA metabolic processes (GO:0006259) and DNA integration (GO:0015074). Membrane (GO:0016020, proteolysis (GO:0006508), and plant-type secondary cell wall biogenesis (GO:0009834) were represented by four counts. Three-count terms were fatty acid biosynthetic process (GO:0006633) and glycosyltransferase activity (GO:0016757). Eight of the two-count transcriptomes and 45 of 1-count transcriptomes were upregulated under blue LED-blue vs. white LED light (Table 8). A total of 42 transcriptomes for GO terms were downregulated. Five-count transcriptomes were obtained for DNA integration (GO:0015074) and binding (GO:0005488). The three-count GO terms included the inositol catabolic process (GO:0019310), RNA-templated DNA biosynthetic process (GO:0006278), and translational initiation (GO:0006413). Two counts were classified into four GO terms, and one count was classified into 33 GO terms (Table 9).

Up-and Downregulation of Transcriptomes in A. membranaceus Sprouts Treated with Red LED-Light vs. White LED-Light
Transcriptomes showing upregulation in red LED vs. white LED comparisons appeared in all 60 GO terms. Among these, the highest count (seven) was obtained for DNA integration (GO:0015074) term. Next, the resolution of meiotic recombination intermediates (GO: 0000712) was represented by six counts. Five counts were classified as proteolysis (GO:0006508), and four counts were grouped as GO terms of flavonoid biosynthetic process (GO:0009813) and defense response (GO:0006952). GO terms with three counts were signal transduction (GO: 0007165), nucleic acid metabolic process (GO: 0090304), and xyloglucan metabolic process (GO: 0010411). In addition, two counts were classified as upregulation of 10 GO terms and one count of 42 GO terms (Table 10). In the downregulation transcriptome result of red LED vs. white LED light, DNA integration (GO:0015074) was represented by six counts. The three-count GO terms were nucleic acid metabolic process (GO:0090304) and binding (GO:0005488), whereas two-count transcriptomes were classified into five GO terms, and one-count transcriptomes were classified into 35 GO terms (Table 11).

Discussion
The quality of light influences the activity of multiple biological pathways in plants and can accordingly have a significant effect on morphological phenotypes. To date, however, there appear to have been no reports regarding morphological changes at the cellular level or a comparative analysis of transcripts in aseptically cultivated plants of the medicinal plant A. membranaceus exposed to different colored LED lights. Previous studies on A. membranaceus transcriptomes have tended to focus on plants grown in general classical cultivation environments [16,17]. In the present study, we analyzed and compared the morphological changes and metabolic mechanisms at the transcriptomic level in A. membranaceus sprouts germinated in vitro under aseptic conditions and illuminated with LED lights of three different wavelengths. The transcriptomic data obtained from this comparative analysis will contribute to gaining a better understanding of the molecular mechanisms underlying the response of A. membranaceus to LED light of different colors, and thereby provide a basis for cultivating medicinal plants with enhanced functional properties.
At the cellular level, we found the in vitro-germinated A. membranaceus sprouts that had been exposed to red LED light were characterized by the largest leaf cells, whereas the thickest stem cells were observed in plants cultivated under white LED illumination (Table 1). Previous studies on the effects of different LEDs on Astragalus growth have found that the growth of A. membranaceus plants cultivated under blue LED illumination was superior to that of plants exposed to either white LED or red LED light when visually measured [18], although these authors detected no significant differences among plants treated with these three light sources with respect to the length and number of leaves. In studies on the effects of LED illumination on the growth of other plants, it has been shown that in Myrtus communis, the highest shoot multiplication effect occurred in plants treated with 5 µM BA(Benzyladenine) grown under red LED light [19], whereas, in Salvia miltiorrhiza, Choi et al. (2020) found red LED light to be more effective than blue LEDs in promoting leaf length and width [20]. Among other studies, it has been reported that blue LED light inhibits vegetative growth, such as seedling growth and root formation, whereas green LED light has been established to have negative effects on plant growth and fresh and dry weights [21,22]. However, although many studies have investigated differences in plant growth responses using LED lights, most assessments of differences in growth have been based on unaided visual evaluations, whereas relatively few studies have examined responses at the cellular level. In this study, using scanning electron microscopy, we examined the leaf cell size and stem thickness of in vitro-germinated A. membranaceus sprouts and accordingly demonstrated red and white LED light sources to be the most effective in promoting leaf cell expansion and stem thickening, respectively. Notably, in this context, the findings of previous studies have indicated that the responses of plants to light of different qualities tend to differ among species. Accordingly, we speculate that our observations for A. membranaceus might also apply to other species in the family Leguminosae. In addition, it has also been reported that exposure to white LED light enhances the antioxidant activity of A. membranaceus [18]. Thus, we assume that the thickening of the stems of A. membranaceus grown under white LED light contributes to enhancing stress resistance and antioxidant activity.
Although there have been several transcriptome analyses using A. membranaceus, the present study is the first to analyze the transcriptome of A. membranaceus sprouts germinated in vitro under sterile conditions. Comparative analysis of plants treated with the three LED light sources revealed that the expression levels of upregulated transcripts in plants exposed to white LED were higher than those in plants cultivated under either red or blue LED sources, indicating that the number of downregulated transcripts is also lowest under white LED light illumination. To functionally annotate the DEGs between different light source treatments, we performed GO analysis. For the blue LED vs. red LED comparison, we detected seven enriched GO sub-categories (phosphorylation, DNA biosynthetic process, binding, membrane, DNA integration, and transport). Similarly, seven enriched categories were detected for the blue LED vs. white LED comparison (binding, RNA-templated DNA biosynthetic process, DNA metabolic process, DNA integration, membrane, proteolysis, and plant-type secondary cell wall biogenesis), whereas for the red LED vs. white LED comparison, transcripts were classified into five categories (DNA integration, resolution of meiotic recombination intermediates, proteolysis, flavonoid biosynthetic process, and defense response).
Conversely, with respect to light source-specific downregulated transcripts, we detected enrichment of DNA integration, an integral component of membrane, and carbon utilization for the blue LED vs. red LED comparison; enrichment of DNA integration and binding for the blue LED vs. white LED comparison; and enrichment of DNA integration for the red LED vs. white LED comparison. Among the upregulated transcripts between blue LED-and red LED-illuminated plants, the most highly enriched GO category was binding (GO:0005488). A GO category that was commonly enriched with upregulated transcripts in all three comparisons is DNA integration (GO:0015074), which is known to function in biological processes that involve the integration of DNA segments into other larger DNA molecules, such as chromosomes. Similarly, the DNA integration category was found to be enriched with downregulated transcripts. We accordingly speculate that each of the three LED light sources has a significant effect on plant cell size and morphogenesis by directly influencing biological processes in sterile A. membranaceus sprouts.
Light plays a particularly important role in the success of in vitro plant tissue culture. LEDs of many artificial light sources are important for plant mass propagation systems in which the color of the light source affects plant growth and development. It has been established that the color of LED light influences the morphogenesis, differentiation, and proliferation of plant cells, tissues, and organ cultures and is essential for the production of secondary metabolites [23,24]. Among the diverse spectrum of secondary metabolites, phenolic compounds have been established to have a broad range of biological properties, including antioxidant, anticancer, antibacterial, and antiallergic activities. As phenolic compounds, flavonoids have been found to have antioxidant and anti-inflammatory properties and play roles in the inhibition of cell division and redox regulation of cells [25,26]. In this regard, an enrichment of phenolic components and high antioxidant capacity has been reported in the leaves and roots of S. miltiorrhiza plants exposed to LED light sources [20]. Furthermore, analysis of phenolic contents in the leaves and roots of Rehmannia glutinosa has revealed that total phenolic contents were highest in those plants exposed to blue LED light, although the highest flavonoid contents were detected in plant cultivated under red LED illumination [27,28]. Changes in the activity of plants promoted by exposure to LED light sources are mediated via the activities of associated genes. In this regard, we found that a GO term of particular interest to us, namely, flavonoid biosynthetic process (GO:0009813), was enriched with four genes that were differentially expressed between red LED-and white LED-illuminated plants. Similar enrichment was detected for the defense response term (GO:0006952), thereby indicating a correlation between plant defense and flavonoid biosynthetic mechanisms.
Flavonoids are a notably abundant class of secondary metabolites present in all terrestrial plants, with more than 10,000 different types believed to occur in different plant species. Flavonoids are a class of phenylpropanoids derived from the shikimate and acetate pathways via the activity of cytoplasmic multienzyme complexes anchored in the endoplasmic reticulum [29]. As defense compounds and developmental regulators, it has been revealed that they can play diverse roles in plant-nematode defense mechanisms by acting as defense compounds or signal molecules that have inhibitory effects on nematodes at different life stages [30]. Therefore, it would be desirable to study the relationships between plant disease defense and various stresses as functions of genes specifically expressed in aseptically cultured A. membranaceus sprouts exposed to white LED light. In addition, it is suggested that the light condition for improving the biomass of A. membranaceus sprouts requires the appropriate use of white light and red light sources.

Preparation of In Vitro-Cultured A. membranaceus Sprouts under Artificial Light Source Conditions
The seeds of A. membranaceus used in this study were purchased from KS Jongmyo (Incheon, Republic of Korea). Seed sterilization was performed to obtain aseptic A. membranaceus sprouts. The seeds were initially placed in 70% EtOH and shaken for 1 min, after which they were transferred to 3% NaClO and shaken again for 5 min. Following five washes with sterile distilled water, the sterilized seeds were transferred to culture bottles containing agar-solidified Murashige and Skoog medium and cultured for 6 weeks. For illumination, we used LED lights of different wavelengths, namely, white (continuous spectrum), red (632 nm, 1.58 µmol/m 2 /s), and blue (465 nm, 1.44 µmol/m 2 /s). The wavelengths of these LED lights were measured using a PG200N illuminometer (United Power Research Technology Co., Zhunan Township, Taiwan). The photoperiod to which the plants were exposed was set to 16 h light and 8 h dark, which was maintained for the 6-week growth period ( Figure 6).

Scanning Electron Microscope Analysis
Leaf and stem samples collected from A. membranaceus plants were fixed with 2% glutaraldehyde and 2% paraformaldehyde in 50 mM cacodylate buffer (pH 7.4) for 1 h at 4 • C [31]. Thereafter, the fixed samples were dehydrated in a graded ethanol series for 10 min and then immersed in mixtures of 100% ethanol and isoamyl acetate (2:1, 1:1, and 1:2), each for 10 min. After immersion in pure isoamyl acetate for 15 min, the isoamyl acetate was removed, and the samples were dried using a critical point dryer. Dried samples were then sputter-coated with a thin layer of gold. Observations were performed at the Korea Basic Science Institute, Chuncheon, using a SUPRA 55VP scanning electron microscope (Carl Zeiss, Oberkochen, Germany) operating at an acceleration voltage of 3 kV. The leaf samples were observed by taking the fourth leaf from the top and measuring 3 leaves.

Scanning Electron Microscope Analysis
Leaf and stem samples collected from A. membranaceus plants were fixed with 2% glutaraldehyde and 2% paraformaldehyde in 50 mM cacodylate buffer (pH 7.4) for 1 h at 4 °C [31]. Thereafter, the fixed samples were dehydrated in a graded ethanol series for 10 min and then immersed in mixtures of 100% ethanol and isoamyl acetate (2:1, 1:1, and 1:2), each for 10 min. After immersion in pure isoamyl acetate for 15 min, the isoamyl acetate was removed, and the samples were dried using a critical point dryer. Dried samples were then sputter-coated with a thin layer of gold. Observations were performed at the Korea Basic Science Institute, Chuncheon, using a SUPRA 55VP scanning electron microscope (Carl Zeiss, Oberkochen, Germany) operating at an acceleration voltage of 3 kV. The leaf samples were observed by taking the fourth leaf from the top and measuring 3 leaves.

RNA-seq Library Construction and Sequencing
Total RNA was extracted from in vitro-cultured A. membranaceus sprouts using Trizol reagent (Invitrogen Scientific, Inc., Waltham, MA, USA), with the purity of the extracted RNA being determined using a microvolume spectrophotometer (Keen Innovative Solutions, Daejeon, Republic of Korea). RNA-seq libraries were constructed using a TruSeq RNA kit (Illumina Inc., San Diego, CA, USA) and sequenced using the Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA, USA).

RNA-Seq Library Construction and Sequencing
Total RNA was extracted from in vitro-cultured A. membranaceus sprouts using Trizol reagent (Invitrogen Scientific, Inc., Waltham, MA, USA), with the purity of the extracted RNA being determined using a microvolume spectrophotometer (Keen Innovative Solutions, Daejeon, Republic of Korea). RNA-seq libraries were constructed using a TruSeq RNA kit (Illumina Inc., San Diego, CA, USA) and sequenced using the Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA, USA).

qPCR Analysis of Reference Genes
cDNA was synthesized using PrimeScript™ RT Master Mix (Perfect Real Time) (Takara Korea Biomedical Inc., Seoul, Republic of Korea) with total RNA used for transcriptome analysis. The qPCR reaction was performed with a CronoSTAR™ 96 Real-Time PCR System (Takara Korea Biomedical Inc., Seoul, Republic of Korea) using TOPreal™ SYBR Green qPCR PreMIX (Enzynomics Co., Ltd., Daejeon, Republic of Korea). qPCR conditions were performed by initial denaturation at 95 • C for 10 min with a volume of 25 µL, followed by Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The data presented in this study are contained within the article.